The effective continuum threshold in dispersive sum rules 
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We study the accuracy of the bound-state parameters obtained with the method of dispersive sum 
rules, one of the most popular theoretical approaches in nonperturbative QCD and hadron physics. 
We make use of a quantum-mechanical potential model since it provides the only possibility to 
probe the reliability and the accuracy of this method: one obtains the bound-state parameters 
from sum rules and compares these results with the exact values calculated from the Schrodinger 
equation. We investigate various possibilities to fix the crucial ingredient of the method of sum 
rules — the effective continuum threshold — and propose modifications which lead to a remarkable 
improvement of the accuracy of the extracted ground-state parameters compared to the standard 
procedures adopted in the method. Although the rigorous control of systematic uncertainties in 
the method of sum rules remains unfeasible, the application of the proposed procedures in QCD 
promises a considerable increase of the actual accuracy of the extracted hadron parameters. 
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1. INTRODUCTION 

The method of dispersive sum rules for the extraction 
of ground-state parameters in QCD was formulated in 
P, Q and since then has been extensively applied to the 
analysis of hadron properties [s] . 

A sum-rule calculation of hadron parameters [11, 0| in- 
volves two steps: (i) one calculates the operator product 
expansion (OPE) for a relevant correlator and formulates 
the sum rule which relates this OPE to the sum over 
hadronic states, and (ii) one attempts to extract ground- 
state parameters by a numerical procedure. Each of these 
steps leads to uncertainties in the final result. 

The first step lies fully within QCD and allows a rigor- 
ous treatment of the uncertainties: the correlator in QCD 
is not known precisely (because of uncertainties in quark 
masses, condensates, as, radiative corrections, etc.) but 
the corresponding errors in the correlator may be con- 
trolled, at least in principle. We refer to such errors as 
the OPE uncertainties. 

The second step is more cumbersome: even if several 
terms of the OPE for the correlator were known pre- 
cisely, the numerical procedures of sum rules should pro- 
vide the range of values which contains the true value of 
the hadron parameter. We call this range the intrinsic 
sum-rule uncertainty. 

In spite of the extensive applications of sum-rules in 
particle physics, including also flavor physics [4], where 
a rigorous error analysis is mandatory, a proper investi- 
gation of the systematic uncertainties of the method has 
been started only recently ^, ifi, 7] . 

The method of sum rules contains a set of prescriptions 
which are believed to allow the control of the accuracy of 
the extracted bound-state parameters (see e.g. Ref. [31). 
The outcome of these prescriptions is claimed to be the 
estimate of the intrinsic sum-rule uncertainty. 



Obviously, the only possibility to acquire an unbiased 
judgement of the reliability of the error estimates in sum 
rules is to apply the method to a problem where the pa- 
rameters of the theory may be fixed and the correspond- 
ing parameters of the ground state may be calculated 
independently and exactly. 

Presently, only quantum-mechanical potential models 
provide such a possibility. A simple harmonic-oscillator 
(HO) potential model, used as a testing ground in 
possesses the essential features of QCD — con- 
finement and asymptotic freedom Q — and has the 
following advantages: (i) the bound-state parameters 
(masses, wave functions, form factors) are known pre- 
cisely; (ii) direct analogues of the QCD correlators may 
be calculated exactly. 

Applying the standard sum-rule machinery, we have 
determined the ground-state decay constant and the 
form factor 6] from the relevant correlators, and con- 
fronted the obtained results with the known exact val- 
ues, probing in this way the accuracy of the method. We 
have clearly demonstrated that the standard procedures 
adopted in the method of sum rules do not yield realistic 
error estimates for the extracted ground-state parame- 
ters. Moreover, we have shown that the uncontrolled 
systematic errors of the form factors are typically much 
larger than those for the decay constants. 

The natural questions which then arise are: (i) Can 
the "standard" procedures of the method of sum rules 
be modified, leading to an improvement of the extracted 
ground-state parameters? (ii) Can one formulate a pro- 
cedure which would provide the interval surely containing 
the actual bound-state parameter? This would mean a 
rigorous control of the intrinsic sum-rule uncertainty. 

In this Letter, we will show that the answer to the 
first question is "yes" , whereas the answer to the second 
question is "no" . 
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The crucial ingredient of sum rules is the effective con- 
tinuum threshold Zcj which governs the accuracy of the 
quark- hadron duality hypothesis, the basic concept of the 
method. We study possible modifications of the standard 
procedure of fixing Zc- In the HO model, relaxing the 
standard assumption of a Borel-parameter independent 
Zc is shown to lead to a significant improvement of the 
extraction of the bound-state parameters, particularly, of 
the form factor. Even though the rigorous control over 
the systematic uncertainties of the ground-state parame- 
ters obtained from sum rules is not feasible (and cannot 
be obtained in principle in problems where the truncated 
OPE is the only input), the application of our findings in 
QCD promises a considerable improvement of the actual 
accuracy of the method. 



2. HARMONIC-OSCILLATOR MODEL 

We consider a non-relativistic HO model defined by 
the Hamiltonian (r = \r\) 

H ^Ho + V{r) , Hq = p^/2m , V{r) = mu;^r^/2 , (1) 

where all features of the bound states are calculable. For 
instance, for the ground (g) state one finds 



i?g = , i?g ^ |vl/g(r = 0)|2 = {mu;/7rf' , 
i^g(g) = exp(-gV4ma;) , (2) 

where the elastic form factor of the ground state is de- 
fined according to (q = \q\) 

F,iq) = i^slJm'f,) = / d^k^ik) ^k-q) , (3) 



with the current operator J(q) given by the kernel 

{r'\J{q)\r) = exp{iq- f) 5^^\f - r') . (4) 

3. POLARIZATION OPERATOR 

In the method of dispersive sum rules the basic quan- 
tity needed for the extraction of the decay constant (i.e., 
of the ground-state wave function at the origin) is the 
correlator of two currents [H. Its quantum- mechanical 
analogue is 

n(r) = (f^ = o|e-^^H = o) , (5) 

where T is the Euclidean time. In the case of the HO 
potential the correlator H(T) is exactly known: 

3/2 

(6) 
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npower(T) = n(r) 
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4. VERTEX FUNCTION 

The basic quantity for the extraction of the form factor 
in the method of dispersive sum rules is the correlator 
of three currents 0]. The analogue of this quantity in 
quantum mechanics is @ 



r(T2,Tl,9) 



{rf^O\e-"^\Jiq)e-"^'\n^O) , 



(7) 



with the operator J{q) being defined in In the HO 
model the exact analytic expression for T(t2,ti, q) was 
obtained in Ref. At equal times ti = T2 = it 
takes the following form: 



r(T,,) ^ H(r) exp(-^tanh(f 



ro(r,q) = Ho(r) exp ( 

' 8m 



(8) 



r(T,q) = r(r,g)-ro(T,q) 



27rT 
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q^uj^ 
24to ' 



In this work we will take into account all the terms in 
the square brackets for both Hpowcr (T) and Fpowor {T, q) . 
Notice that each term is a power in T and/or q^. Thus, 
retaining a fixed number of power corrections restricts the 
convergence of Hpower(I^) and T power {T,q) to the region 
of not too large values of T and/or q"^, as it happens in 
QCD when the OPE series is truncated. 



5. GROUND-STATE PARAMETERS 

Making use of the quark-hadron duality hypothesis, 
which assumes that the excited-state contribution is dual 
to the high-energy region of the free-quark diagrams, one 
gets the sum rules for Rg 



i?ge-^«^ =Hpowor(T)-f J dzpo{z)e-^^, (9) 



and for the form factor Fg(q) 



(10) 



RgFg{q)e-^-' =rpower(r,(z) 

+ j dzi j dz2e--^^''+^^^^Ao{zi,Z2,q) , 



where po{z) and Ao(2:i, Z2, q) are the known spectral den- 
sities of the two- and three-point Feynman diagrams of 
the non-relativistic field theory [3, [9| . 

The relations (|9ll0p constitute the definitions of the ex- 
act effective continuum thresholds z^(T) and Zctf{T,q). 
Their full T- and g-dependences can be obtained by solv- 
ing Eqs. (|9ll0p using the exact bound-state parameters 
-Rg and Fg{q) as well as the exact power expansions 
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npowor(T') and Tpowcr(T, q). In the HO model this can be 
easily done numerically. Without loss of generality we set 
m = Lo and show the corresponding results in Fig. [TJ It 
can clearly be seen that the effective continuum threshold 
Zcs{T, q) does depend upon both T and q. 
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Fig. 1; Effective continuum threshold Zcs{T, q) for the 3-point 
function, obtained by solving numerically Eq. (|10|l using the 
exact bound-state parameters Tig and Fg{q) as well as the 
exact power expansion rpower(r, q), versus the Euclidean time 
T at fixed values of the momentum transfer q (a) and versus 
q at fixed values of T (b). In (a) the vertical dashed lines 
identify the fiducial range in T (see text), while the solid lines 
are linear fits of Zcff (T, g) in the fiducial range. 

Let us consider a restricted problem when the energy 
Eg of the ground state is known, and try to determine 
its elastic form factor from the sum rule pl}|) . 

First, according to [l| we should determine the Borel 
window (or the fiducial range), where the sum rule may 
be used for the extraction of the ground-state parameter: 
i) the lower boundary of the T-window is found from the 
requirement that the ground state gives a sizable (we 
require more than 50%) contribution to the correlator; 
and ii) the upper boundary of the T-window is obtained 
from the condition that the truncated OPE gives a good 
approximation to the exact correlator. Since in the HO 
model the power corrections are exactly known, the up- 
per boundary is T = oo. However, to be close to realistic 
situations when only a limited number of power correc- 
tions is available, we take from our study of Ref. 01 the 



fiducial range 0.7 < wT < 1.2 (see Fig.ITTa)) [l^. 

Second, we must choose a criterion to approximate the 
effective continuum threshold Zcs{T,q). In this work we 
compare three different approximations: 

Z,siT,q) « z^{q) , (11) 
Zeff(T,q) « z^iq) + z^iq) cjT , (12) 
Zcff(T,<7) « z^{q) + z?iq)LuT + z§iq)LO^T^ .{13) 

The standard procedure adopted in the sum-rule method 
is to assume a T- independent value, i.e. Eq. (Ill|l . 

At each value of q we fix the parameters appearing on 
the r.h.s of Eqs. (fTTHT5|) in the following way: we define 
the dual energy, £duai(r, g), as 



EduediT,q) 



dT 



logrdual(T,q,Zcff(r,g)) , (14) 



where Fduai is the r.h.s. of Eq. (|T0|) calculated using the 
approximations (|llH13p for Zoff(T, q). Then we calcu- 
late Eduai{T, q) at several values oi T = Ti {i = 1, . . . , N) 
chosen uniformly in the fiducial range and finally we min- 
imize the squared difference with the exact value E„: 



X 



1 ^ 



{T,,q)-Eg]' 



(15) 



The results for the dual form factor i^duaK*?), obtained 
via Eqs. @ and PU)) using the correlator Fduai after op- 
timizing the parameters of the three approximations ([Tl]~ 
[T5)) . are shown in Fig. [21 Note that because of current 
conservation the form factor should obey the absolute 
normalization Fduai (g = 0) = 1. We therefore require 
z^g{T) = z^s{T,q = 0), so that the r.h.s. of Eqs. (PTU]) 
coincide and the dual form factor is properly normalized. 

1.2 



1.1 



1.0 



0.9 



• 


constant 




□ 


linear 


□ 


♦ 


quadratic 


n 



0.8 

0.0 0.5 1.0 1.5 2.0 2.5 3.0 
q / CO 

Fig. 2: Ratio of the dual form factor F^uaiiq), extracted 
from the sum rule (|10p using different approximations for 
Zcs{T, q), and the exact ground-state form factor i^g(g), given 
by Eq. The dots, squares and diamonds correspond, re- 
spectively, to the results obtained using the constant 
the linear p2p and the quadratic (|13|) approximations for the 
T-dependence of Zctt (T, q) . 

Let us consider first the case of the T-independent ap- 
proximation pip . The criterion of minimizing the dlB 
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leads to ZQ{q) for which the dual energy E^uai{T,q) dif- 
fers from Eg by less than 0.1% and to the dual form fac- 
tor Fduediq) which is largely T-independent in the whole 
fiducial range. Such a stability, usually referred to as the 
Borel stability, is often (erroneously) claimed to be the 
way to control the accuracy of the extracted form fac- 
tor. From Fig. [2] it can be seen that the T-independent 
approximation (fTT|) works well (better than 2%) at low 
values of q. However, for q > 1.7 lo the T-independent 
ansatz does not work at all, since it is impossible to re- 
produce the ground-state energy in the fiducial range. 

When the linear approximation for Zcs{T,q) is 
considered, the form factor can be extracted also for q > 
1.7 Lu. Note that the exact effective continuum threshold 
can be very well approximated by a linear function of T 
in the whole fiducial range [see Fig. [IJa)]. Nevertheless, 
this is not a guarantee that one can extract the exact form 
factor: deviations of the order of several percent can be 
produced after minimization of Eq. (jlSp up to q « 2uj 
and uncertainties of the order of 10 20% may plague 
the extracted form factor at g > 2a;. 

One may try to go further and consider the quadratic 
Ansatz for ZceiT,q) (jl3p . However, as can be seen from 
Fig. [21 this leads to certain instabilities in the extracted 
value of the form factor. These instabilities just reflect 
the fact that the unique solution to the problem of ex- 
tracting the form factor from the correlator in a lim- 
ited T-window does not exist Therefore, there is no 
way to get a systematic improvement in the accuracy of 
the extracted form factor by increasing the degree of the 
polynomial Ansatz for Zeg{T^ q). 

Nevertheless, it is worth emphasizing that in the HO 
model the comparison between the form factor extracted 
assuming Eq. and the one obtained using Eq. 
gives a realistic estimate of the accuracy in a wide range 
of values of q. Whether this feature persists in QCD is 
an interesting and important issue to be addressed in the 
future. 



6. CONCLUSIONS 

Let us summarize the main messages of our analysis: 
• The knowledge of the correlator in a limited range of 



relatively small Euclidean times T (that is, large Borel 
masses) is not sufficient for the determination of the 
ground-state parameters. In addition to the OPE for the 
relevant correlator, one needs an independent criterion 
for fixing the effective continuum threshold. 

• Assuming a T-independent (i.e., a Borel-parameter 
independent) effective continuum threshold the error of 
the extracted ground-state parameter (both decay con- 
stant and form factor) turns out to be typically much 
larger than (i) the error of the description of the exact 
correlator by the truncated OPE and (ii) the variation 
of the bound-state parameter in the fiducial range (i.e., 
Borel window) . The latter point is of particular relevance 
since the Borel stability is usually believed to control the 
accuracy of the extracted ground-state parameter. Ob- 
viously, this is not the case (see also Refs. d, (EH)- 

• Allowing for a T-dependent effective continuum 
threshold and fixing it according to Eq. leads to 
evident improvements in the extracted ground-state pa- 
rameters. This was shown for the decay constant and the 
form factor in the HO model. Moreover, in this model 
the variation of the form factor extracted using differ- 
ent approximations for z^g gives de facto a realistic error 
estimate. Unfortunately, the use of higher polynomial 
approximations leads to instabilities in the fitting proce- 
dures. It is, therefore, impossible to construct a system- 
atic procedure which would converge to the exact effec- 
tive continuum threshold. As the result, rigorous error 
estimates cannot be obtained. 

The impossibility to get a rigorous control over the sys- 
tematic errors of the extracted ground-state parameters 
is the weak feature of the method of sum rules and an 
obstacle for using the results from sum rules in problems 
where rigorous error estimates are required. 

In spite of this weakness, the application of the pro- 
posed modifications of the method in QCD seems very 
promising and may lead to a considerable increase of 
the actual accuracy of the calculated hadron parameters. 
This issue deserves a serious investigation. 
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